Spontaneous Four- Wave Mixing of de Broglie Waves: Beyond Optics 
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We investigate the atom-optical analog of degenerate four-wave mixing of photons by colliding 
two Bose-Einstein condensates (BECs) of metastable helium and measuring the resulting momentum 
distribution of the scattered atoms with a time and space resolved detector. For the case of photons, 
phase matching conditions completely define the final state of the system, and in the case of two 
colliding BECs, simple analogy implies a spherical momentum distribution of scattered atoms. We 
find however, that the final momenta of the scattered atoms instead lie on an ellipsoid whose radii 
are smaller than the initial collision momentum. Numerical and analytical calculations agree with 
the measurements, and reveal the interplay between many-body effects, mean- field interaction, and 
the anisotropy of the source condensate. 
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The field of atom optics has developed to the point that 
one can now speak of the beginning of "quantum atom 
optics" [l[ in which atoms are manipulated in ways simi- 
lar to photons and in which quantum fluctuations and en- 
tanglement play an important role. The demonstration of 
atom pair production 0, Q , either from the dissociation 
of ultra-cold molecules ,a process analogous to paramet- 
ric down-conversion [JQ > or from collisions of BECs 0- 
[Io| , analogous to four- wave mixing (FWM) [lll-[21|, holds 
considerable promise for generating atomic squeezed 
states and demonstrating nonlocal Einstein-Podolsky- 
Rosen (EPR) correlations @, H, IH, H3| . In both these 
systems, atom-atom interactions play the role of the non- 
linear medium that allows conversion processes. Atoms 
are not, however, exactly like photons, and in spite of 
their formal similarity, the processes of pair production 
of photons and of atoms exhibit some interesting and 
even surprising differences that must be understood in 
order for the quantum atom optics field to advance. In 
this work, we discuss one such effect. 

In optical FWM or parametric down conversion {24j|, 
energy conservation requires that the sum of the energies 
of the outgoing photons be fixed by the energy of the 
input photon(s). Phase matching requirements impose 
constraints on the directions and values of the individual 
photon momenta. A simple case is degenerate, sponta- 
neous FWM (i.e. two input photons of equal energy) in 
an isotropic medium, for which energy conservation and 
phase matching require that the momenta of the output 
photons lie on a spherical shell whose radius is that of 
the momenta of the input photons. 

We have performed the atom optical analog of degen- 



erate FWM in colliding BECs while paying careful at- 
tention to the momenta of the outgoing atoms. We find 
that unlike the optical case, the output momenta do not 
lie on a sphere, but rather on an ellipsoid with short ra- 
dius smaller than the input momentum. This behavior is 
due to a subtle combination of atom- atom interactions, 
which impose an energy cost for pair production, and the 
anisotropy of the condensates, which affects the scattered 
atoms as they leave the interaction region. 

Although an analogous effect could exist in optics, op- 
tical nonlinearities are typically so small that the effect is 
negligible. However in the process of high-harmonic gen- 
eration in intense laser fields, a similar effect has been 
discussed [H[. There, phase matching conditions can 
become significantly intensity dependent, and the pon- 
deromotive acceleration of electrons alters the phase and 
energy balance of the harmonic generation process. Thus 
the ponderomotive force plays a role loosely analogous to 
that of the mean-field repulsion in our problem. 

To fully understand the results, we have simulated the 
BEC collision using a fully quantum, first-principles nu- 
merical calculation based on the positive- P representa- 
tion method 1?], EqI, and find quantitative agreement 



with the experiment. We have also analyzed the prob- 
lem using a stochastic implementation of the Bogoliubov 
approach, which allows us to identify and illustrate the 
contributions of various interaction effects in the process. 

The experimental setup is similar to that described 
in jjjj. We start from a BEC of ~ 10 5 atoms magnetically 
trapped in the m x = 1 sublevel of the 2 3 Si metastable 
state of helium-4. The trap is cylindrically symmetric 
with axial and radial frequencies of 47 Hz and 1150 Hz, 
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FIG. 1: (Color online) (a) Geometry of the Bragg beams and 
level scheme of the 2 3 Si - 2 3 P transition of 4 He (at 1083 
nm). A Bragg pulse of two 7r-polarized laser beams (shown 
by the two arrows) detuned by A/2tv = 600 MHz produces 
two counterpropagating BECs that separate along their ra- 
dial dimension at approximately 45° to the vertical (z) axis 
at relative velocity 2vo. (b) Schematic diagram of the colli- 
sion geometry in the center-of-mass frame in which we denote 
the collision axis as Z. The two disks represent the colliding 
condensates in momentum space. The sphere represents the 
halo of scattered atoms. The cigar shaped initial condensate 
with axial direction X = x is shown in the center. We analyze 
the experimental data in the Xy-plane. 



respectively. The bias field of ~ 0.25 G along the x-axis 
defines the quantization axis. 

To generate the two colliding BECs, we use a two-step 
process. First, the atoms are transferred to the m x = 
state by a stimulated Raman transition. Using a 4 /is 
long pulse, we transfer 90% of the atoms to this magnet- 
ically untrapped state. 1 /is after the end of the Raman 
pulse, the BEC is split into two counterpropagating con- 
densates with a Bragg pulse driven by two laser beams 
propagating at approximately 90°, as shown in Fig. [I] (a). 
The parameters of the Bragg pulse are adjusted to trans- 
fer half of the atoms to a state moving at relative velocity 
2vo in the ^/z-plane, with = 7.31 cm/s, which is ~ 4 
times the speed of sound in the center of the BEC. The 
condensates thus separate along the radial axis, unlike 
in the experiment of Ref. Q. To analyze the data we 
will use a center-of-mass reference frame, in which the 
collision axis is defined as Z (tilted by about 45° from 
z), X = x, and Y is orthogonal to Z and X (see Fig. [1]). 

After the collision, the atoms fall onto a microchan- 
nel plate detector placed 46.5 cm below the trap center. 
A delay line anode permits reconstruction of a 3D im- 
age of the cloud of atoms. The flight time to the detec- 
tor (300 ms), is long enough that the 3D reconstruction 
gives a 3D image of the velocity distribution after the 
collision. Binary, s-wave collisions between atoms in the 
BECs should (naively) result in the scattered particles 
being uniformly distributed on a sphere in velocity space 
with radius equal to the collision velocity vq. The col- 
lision along the radial axis allows access to the entire 
collision halo in a plane containing the anisotropy of the 
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FIG. 2: (Color online) (a) Average momentum space density 
n(kx,ky) (in arb. units, from ~ 1500 experimental runs) 
of the experimentally observed scattering halo on the equa- 
torial plane (kx->ky); the density is averaged over a disk of 
thickness [— O.lfco, +0.1fco] along kz- (b) Same as in (a) but 
from the positive- P simulation (see text) after 70 /is collision 
time, in units of 10 -18 m 3 . (c) Plot of the peak radius of the 
scattering halo on the equatorial plane versus the azimuthal 
angle cj). Black squares are experimental data, while the red 
circles are from the simulation. The data is binned into 18 
angular bins of Acj) — 20°, and each data point for the peak 
radius is derived from a Gaussian fit to the radial distribution 
n(kR,<p) = n(kx — kR cos 0, ky — kRsm(j>) at the respective 
angle <j) (the error bars show the statistical uncertainty in the 
fits; in addition, there is a systematic uncertainty of ±1.5% in 
the determination of the average radius of the sphere). The 
smooth line is a sinusoidal fit to the experimental data. 

BEC (the XF-plane) without distortion from the con- 
densates. As in Ref. Q, we observe a strong correlation 
between atoms with opposite velocities confirming that 
the observed halo is indeed the result of binary collisions. 

In Fig. [2] (a) we show a slice of the scattering halo in the 
XF-plane that reveals its annular structure. A dashed 
circle of radius 1, indicating the momentum hko = mvo, 
is shown for comparison. We can see that the ring cor- 
responding to the mean momentum of scattered atoms 
does not lie exactly on the dashed line, but rather slightly 
within it, and that the deviation is anisotropic. The ring 
thickness and density are also anisotropic, though in the 
present work we concentrate on the behavior of the ra- 
dius. To analyze the data more quantitatively, we divide 
the ring into azimuthal sectors and fit a Gaussian peak 
plus a linearly sloped background to extract a value for 
the halo radius as a function of the angle <\> [20]. It is 
clear from Fig. [2] (c) that the radius of the halo in mo- 
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mentum space varies approximately sinusoidally by ±2% 
and that it is almost always smaller than k$. 

To understand this result qualitatively, we first con- 
sider the energy balance for pair production in a homo- 
geneous BEC. Removing an atom from the condensate 
liberates an energy corresponding to the chemical poten- 
tial, gp, where g = Anah 2 /m, a is the s- wave scatter- 
ing length, and p the density. Here, we have two coun- 
terpropagating condensates (each having density p/2), 
which for simplicity we model as plane waves. In the 
presence of the spatial modulation due to their interfer- 
ence, the energy liberated by removing one atom changes 
to 3gp/2 [26]. On the other hand, placing an atom in a 
scattering mode requires an energy 2gp since the scat- 
tered atom is distinguishable from those in the conden- 
sate. Energy conservation, including the mean-field con- 
tributions, gives 



first-princi ple s positive- P simulations similar to those in 
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where we denote the absolute momentum of one scat- 
tered atom hk s . Thus, the initial scattered momentum is 
smaller than the ingoing momentum, k s < ko. This ef- 
fect was observed in a numerical simulation in Ref. [l4| ; a 
similar effect was discussed in Ref. [s[ . Using plane waves 
to model the BECs is of course a crude approximation, 
but if we replace p by the central density of an inhomo- 
geneous BEC, we find k s = 0.96 ko for the experimental 
parameters. 

In addition to this initial energy balance analysis, a 
second effect must be taken into account. Once created, 
the scattered atoms escape from the condensate region 
and gain energy from the mean-field interaction poten- 
tial. The effect is similar to that reported in Ref. [26| , an 
experiment which observed the mutual repulsion of two 
BECs after Bragg diffraction. If the source BEC were sta- 
tionary, atoms would gain a kinetic energy 2gp as they 
roll-off the mean-field potential. In our system however, 
the potential also evolves in time and goes to zero in the 
XF-plane on a timescale corresponding to the time for 
the two condensates to separate (~ 70 ps). The rapid 
vanishing of the potential on the equatorial plane has a 
very different effect on scattered atoms moving in the X 
and Y directions. Atoms moving along Y, the small di- 
mension of the trap, escape the condensate overlap region 
on a timescale of ~ 40 /is, faster than the condensates 
can separate. As a result, these atoms are accelerated by 
a steep potential gradient and regain part of the energy 
2gp (part - because the potential itself is reduced dur- 
ing the separation). On the other hand, atoms moving 
along X, the long axis of the trap, do not escape before 
the condensates separate and thus experience much less 
acceleration. Accordingly the observed momentum along 
the X direction is smaller than along Y", and much closer 
to the shifted value predicted by Eq. (pp). 

To describe this experiment quantitatively we perform 



Refs. [17|, |20|. Here, the multimode dynamics of the 



atomic field operators \I/(x, t) and \J^(x, ^ f or the m x = 
state is fully modeled by two independent complex c- 
fields, \I/(x, t) and \I/(x, £), satisfying the Ito stochastic 
differential equations: 



ihd t ^(^t) = ^p(*,*)* + v*ft0*Ci(x,*), (2) 
-ifift*(x,t) = A}p(^)^ + \/^%^C2(x,£). 

Here, ^L G p(^,^) = -ft 2 V 2 /(2ra) + g^ is a deter- 
ministic part similar to the mean-field Gross-Pitaevskii 
(GP) equation, Cj( x ? t) (j — 1, 2) are real indepen- 
dent noise sources with zero mean and correlations 
(0(x,£)a(x',t')> = S jk 5W(x - x')5(t - t'), while g = 
4:7rh 2 a/m uses a = 5.3 nm [3] for the rn x = atoms. 

The initial condition for the outcoupled BEC in the 
m x = state (assuming perfect outcoupling for sim- 
plicity) is a coherent state with the same density pro- 
file p(x) as the trapped BEC in the m x = 1 state, with 
a = 7.51 nm [13, A^o = 10 5 atoms. Modulating this with 
a standing wave imparts initial momenta ±fco in the Z 
direction, 



(x, 0) = (*(x, 0)) = VpWT^ {e ikoZ + e- ikoZ ) 



(3) 



and models the Bragg pulse that splits the BEC into two 
equal halves described in the center-of-mass frame. The 
initial density p(x) is obtained as the ground state so- 
lution to the GP equation in the trap, and \I/(x, 0) = 
\I/(x, 0)*. The results of this simulation are shown in 
Fig. [2] (b) and (c) for t = 70 ps at which time the con- 
densates have fully separated and the collision is over. 
The result of the simulation is in reasonable agreement 
with the experiment. The remaining discrepancy could 
be because the experiment, unlike the simulation, aver- 
ages over a broad distribution of initial atom numbers. 
Since large condensates scatter more atoms, these events 
have more statistical weight and bias the data towards 
larger modulations. 

In order to confirm the qualitative mean-field mech- 
anisms described above, we also perform an analysis 
of the collision dynamics using a time-adaptive Bogoli- 
ubov approach [28|, in which the atomic field operator 
is split into the mean-field (^o) and fluctuating compo- 
nents, $(x,t) = ^o(x,£) +£(x,t). The coherent BEC 
wavefunction ^o(x, t) evolves according to the standard 
time-dependent GP equation, with the initial condition 
given by Eq. (j3j). The fluctuating component <5(x, t) de- 
scribes incoherent scattered atoms, and is initially in the 
vacuum state. In the Bogoliubov approach, 5 evolves as 



ihdtSfa *) = Wo(x, t)S + 0(x, t)5l 



(4) 



Here, H (x,t) = -h 2 V 2 /{2m) + 2^|^o(x, t)\ 2 contains 
the kinetic energy and the mean-field potential energy 
2g|^o(x, t)\ 2 for scattered atoms. The effective coupling 
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FIG. 3: (Color online) Predictions for the peak radius of 
the scattering halo as in Fig. [2] (c), after the end of the 
collision (72 fis), with various controlled changes. Red-*: 
full positive-P calculation, Eq. (|2} [same as in Fig. [2] (c)]; 
Black-E\: anisotropic Bogoliubov calculation, Eq. J5}; Blue- 
<: anisotropic Bogoliubov, but with mean- field potentials 
oc g\ipo\ 2 removed from Eq. (j5} and from the GP equation for 
^>o(x, t); Green-x: full Bogoliubov, but with spherical BECs 
and unchanged peak density p(0) (200 /jls). 



context of dissociation of a BEC of molecular dimers, 
which produces atom-atom correlations similar to four- 
wave mixing. In addition to the kinematic effects we re- 
port here, mean-field effects will also affect the phases of 
the associated two-particle wavefunctions. Future work 
must carefully evaluate the effects of such (anisotropic 
and possibly fluctuating) phase shifts on observables like 
the contrast of one- and two-particle interference fringes. 
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£/(x, t) = g^o(x, t) 2 causes spontaneous pair production 
of scattered atoms. The dynamics of the field S is then 
formulated using the positive-P representation [28], lead- 
ing to the (stochastic field) evolution equations 



ihd t S(x, t) 
-ihdtSfa t) 



n s- 
n s- 



gs + v#Ci(x,t), 



(5) 



which, unlike the full calculation (j2j), are stable in time 
because the noise is non- multiplicative. This method 
takes into account the temporal evolution and spatial 
separation of the two condensates; the stochastic for- 
mulation of the evolution of the field <5(x, t) makes ex- 
plicit diagonalizations on the (enormous) Hilbert space 
unnecessary. As condensate depletion is ~ 1.5% here, 
the stochastic Bogoliubov results are in excellent agree- 
ment with the positive-P simulations, as seen in Fig. 03 

Figure [3] also shows simulations performed with con- 
trolled changes applied to the system. The green (x) 
points use a spherical initial condensate and show no 
anisotropy in the scattering sphere, unlike the black (□) 
squares for the anisotropic case. The blue (<) points have 
no mean-field potential, confirming that this potential is 
essential for both the radius shift and the ellipticity. 

The ability to detect three dimensional momentum 
vectors of individual atoms allows the identification 
of small, previously unseen anomalies in the scatter- 
ing "sphere" resulting from a simple collision between 
two condensates. First-principles simulations reproduce 
these small anomalies and help us to identify the im- 
portant physical processes. An important application of 
pair production is the study and exploitation of quantum 
correlations between the pairs, for example via Bell and 
EPR type experiments 29|, [Soj . A matter- wave analogue 
of the optical EPR experiment with parametric down- 
conversion [30| has been discussed in Ref. [22| in the 
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